o 
o 

(N 
Or 



> 
o 

On 
O 

o 
S3 

9 L,: 
6 

CO 
C3 



X 



A Measurement of the UHECR Spectrum with the HiRes FADC Detector 
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We have measured the energy spectrum of ultra-high energy cosmic rays (UHECR) with the HiRes FADC 
detector (HiRes-2) in monocular mode. A detailed Monte Carlo simulation of the detector response to air showers 
has been used to calculate the energy dependent acceptance of the air fluorescence detector. The measured 
spectrum complements the measurement by the HiRes- 1 detector down to lower energies. Systematic effects of 
the assumed input spectrum and composition on the aperture are presented, as well as systematics due to the 
atmosphere. 



1. Introduction: Mono vs. Stereo 

The two air fluorescence detectors of the HiRes 
experiment provide stereoscopic observation of 
extensive air showers (EAS), initiated in the at- 
mosphere by ultra-high energy cosmic rays. The 
stereo view leads to an improved resolution in the 
reconstructed arrival direction and energy of the 
primary cosmic ray particle due to better con- 
straints on the location of the shower axis. How- 
ever, even if observation with both "eyes" is the 
preferred method, there are good reasons for the 
analysis of monocular data. Since the HiRes- 1 de- 
tector began operation two years before HiRes-2, 
the HiRes- 1 dataset provides larger statistics than 
the stereo data. This is especially significant in 
the measurement of the cosmic ray energy spec- 
trum at energies around the GZK flux suppres- 
sion. 

The HiRes-2 monocular dataset, which will 
be presented here, offers a different advantage: 
It contains well reconstructed events at energies 
lower than HiRes-1 mono and stereo events. The 
HiRes-2 detector uses FADC electronics to record 
signals at a frequency of 10 MHz, which leads 
to a better time resolution than is achieved with 
the sample- and-hold electronics of the HiRes-1 
detector. Furthermore, the two rings of mirrors 
of HiRes-2 provide twice the coverage in eleva- 
tion angle that is achieved with the single ring of 



mirrors of the HiRes-1 detector. Events recorded 
with HiRes-2 can be reliably reconstructed down 
to energies of about 10 17 eU. This is also about 
a decade lower in energy than the lower limit for 
stereo events, which is constrained by the separa- 
tion between the two detectors: Events with the 
lowest observable energy lie half-way between the 
two detectors and have thus a distance of about 
6 km from each detector, which sets a lower limit 
to their observable energies. 

The HiRes-2 measurement of the UHECR spec- 
trum is intended to complement the HiRes- 1 mea- 
surement at lower energies, where the features of 
the "second knee" and the "ankle" provide rec- 
ognizable signatures of the cosmic ray spectrum, 
which can be used to compare the HiRes mea- 
surements to other experiments. 

2. Monte Carlo Simulation 

The calculation of the cosmic ray energy spec- 
trum from the measured event distribution is a 
problem of unfolding the true spectrum of cos- 
mic rays at their arrival in the earth's atmosphere 
from the distorsions of the detector response, i.e. 
the acceptance of the detector and its limited res- 
olution. 

The purpose of the Monte Carlo (MC) simula- 
tion in the HiRes experiment is to give an accu- 
rate description of the geometry and energy de- 
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Figure 1. top: Distribution of photoelectrons per 
degree of track for data (filled squares) and MC 
(open squares). Histograms have been normal- 
ized to cover the same area, bottom: Ratio of 
data over MC with linear fit. 



pendent acceptance of the detector and of its res- 
olution. The simulation consists basically of two 
parts: an air shower generator and a detector re- 
sponse MC. In the first part of the simulation, 
large sets of air showers are generated with sev- 
eral discrete energies and with different primary 
particles, using CORSIKA [T| and QGSJet g|. 
Their profiles are saved in a library of air show- 
ers. From this library, the individual profiles can 
be read into the detector response MC and used 
to simulate EAS at different geometries. Ener- 
gies in the detector response MC are chosen from 
a given continuous input spectrum. Profiles from 
library showers generated at a nearby, discrete 
energy are scaled to the chosen energy. 

The detector response simulation includes: 
generation of fluorescence and Cerenkov light at 
the shower, propagation of light through the at- 
mosphere, ray tracing of photons through the op- 
tical path of the detector, PMT response to the 
signal, simulation of noise, electronics and trigger 
simulation. Events accepted by the trigger are 
written out in the same format as the data. 



Two databases have been generated and are 
read by the detector response MC in order to al- 
low event simulation under the exact data tak- 
ing conditions. A trigger database contains in- 
formation on the livetime of the detector, dead 
mirrors and variable trigger settings. An atmo- 
spheric database provides hourly measurements 
of the aerosol content of the atmosphere. An av- 
erage atmosphere has been used in this analysis 
to allow a direct comparison with the HiRes-1 
measurement. 
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Figure 2. Exposure as a function of energy, cal- 
culated from MC simulations. 



3. Event Reconstruction and Analysis 

As the air shower develops in the atmosphere, 
its image sweeps across the photomultiplier tube 
clusters, causing a track of triggered phototubes. 
The first step in the reconstruction of an event 
is the determination of the shower-detector plane 
from a fit to this track. The geometry of the 
shower axis within this plane and the distance 
of the shower from the detector are found in the 
HiRes-2 monocular reconstruction with a fit of 
the times of triggered phototubes versus their an- 
gle along the shower track. Once the location 
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and geometry of the shower axis have been de- 
termined, the charged particle profile of the air 
shower is reconstructed from the distribution of 
recorded FADC counts in triggered phototubes. 
Absolute calibration of the photomultiplier tubes 
with a calibrated light source and measurements 
of the aerosol density of the atmosphere with 
steerable lasers ensure the accurate calculation of 
the shower profile from the FADC signals. Us- 
ing the atmosphere as a calorimeter, we can cal- 
culate the total deposited energy of the shower 
from the charged particle profile after subtraction 
of Cerenkov light. Adding about 10% for "miss- 
ing energy" , which is deposited into the ground 
in muons and neutrinos or lost in nuclear excita- 
tions, yields the total energy of the primary cos- 
mic ray particle. 

The same quality cuts and the same reconstruc- 
tion procedure are applied to both data and sim- 
ulated events, which allows us to directly com- 
pare all the details of our MC simulation against 
the actual experiment. An extensive set of data- 
MC comparison plots — including e.g. the distance 
of air showers, their angular distributions, num- 
ber of triggered mirrors and phototubes, recon- 
structed profiles and energies, etc. — tests how 
well our MC simulates the actual detector. An 
example for the overall very good agreement be- 
tween data and MC is shown in Figure ^ The 
distributions of photoelectrons per degree of track 
are measures of the amount of light that is seen 
in the data and generated in the MC. The simu- 
lation closely follows the data distribution. More 
details on the HiRes-2 analysis can be found in 0. 

With our MC simulation proven to be realistic, 
we can calculate the acceptance of the HiRcs-2 
detector, which is given by the ratio of fully recon- 
structed events (Rmc) over au generated events 
(Gmc) in each energy bin. Multiplying the ac- 
ceptance by the geometrical aperture AVL and the 
livetime t of the detector (about 531 hours of an- 
alyzed data) yields the exposure of the detector. 
We fit the simulated exposure to an appropriately 
chosen function in order to smooth out statistical 
fluctuations. The result is shown in Figure 

The differential flux J in each energy bin can 



now be calculated as: 
J(Ei) = N(Ei) *^=* -5 7— 
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where N(Ei) is the number of data events in the 
energy bin and AE is the binwidth. 
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Figure 3. Energy spectrum (J * E 3 ) measured by 
the HiRes-2 detector in monocular mode. 



4. The HiRes-2 Spectrum 

The energy spectrum measured by the HiRes-2 
detector on clear, moonless nights from Decem- 
ber 1999 to September 2001 (~ 2700 events) is 
shown in Figure |3| We observe the " ankle" at 
about l0 1S 5 eV. Our updated measurement is 
consistent with the GZK flux suppression, but 
statistics at the high energy end are low. Fig- 
ure 0] shows the spectra measured by HiRes-1, 
HiRes-2 and the Fly's Eye stereo experiment su- 
perimposed. The two monocular HiRes spectra 
agree with each other and with our predecessor 
experiment. There is a slight difference in the 
slopes of the HiRes and Fly's Eye spectra above 
the " ankle" , which is most likely due to a dif- 
ferent atmospheric model used in the Fly's Eye 
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experiment. The low energy end of the HiRes-2 
spectrum is consistent with the Fly's Eye spec- 
trum. However, due to low statistics, HiRes can- 
not claim detection of the " second knee" , seen by 
Fly's Eye at 10 176 eU. 
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Figure 4. Energy spectra (J * E 3 ) measured by 
HiRes-1 (squares), HiRes-2 (circles) and Fly's Eye 
stereo (triangles). 



mated by calculating the acceptance for different 
model assumptions in the MC • We have varied 
the shape of the input spectrum, from which MC 
event energies are chosen. Figure [S] shows a plot 
comparing the energy distributions of a MC set 
using a simple E~ 3 power law as input with about 
half of our analyzed data. The disagreement of 
the data with this assumption can be seen clearly 
from the tilted ratio plot. 
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5. Studies of Systematics 

The main systematic uncertainties in the HiRes 
monocular spectra have been reported in |3] : Un- 
certainties in the absolute phototube calibration, 
fluorescence yield, "missing energy" correction 
and atmospheric calibration add up to a total un- 
certainty in the flux of ±31%. Here, we want to 
discuss a potential impact on the calculated ac- 
ceptance from our assumptions about the MC in- 
put energy spectrum and input composition, and 
the assumption of an average atmosphere in our 
simulation and analysis. 

5.1. Input Spectrum 

Using MC simulations to determine the accep- 
tance of the detector could potentially introduce 
a bias into the analysis due to a wrong model 
used in the MC. This potential bias can be esti- 



Figure 5. top: Energy distribution of real 
events (filled squares) and simulated events (open 
squares) using an E~ 3 input spectrum, bottom: 
Ratio of data over MC. 



Our regular MC uses an input spectrum whose 
shape fits the Fly's Eye stereo spectrum below 
the "ankle" and follows a straight line fit to the 
HiRes- 1 data for higher energies. The GZK fea- 
ture has not been included. A comparison plot for 
this input spectrum is shown in Figure El Here, 
the agreement is very good. 

We have calculated acceptances for both MC 
sets. The bias we are avoiding by using a realis- 
tic input spectrum is of the order of ±20%. This 
bias would translate directly into the measured 
spectrum, if one assumed a wrong (i.e. E ) in- 
put spectrum in the MC simulation. 
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Figure 6. top: Energy distribution of real events 
(filled squares) and simulated events using a re- 
alistic input spectrum as described in the text 
(open squares) . bottom: Ratio of data over MC. 



5.2. Input Composition 

We determine the fraction of showers initiated 
by light and heavy, i.e. proton and iron, cosmic 
rays in the MC from composition measurements 
by the HiRes/MIA |H| and HiRes stereo [7] exper- 
iments. 

The uncertainties in the HiRes/MIA measure- 
ment that translate into our spectrum calculation 
add up to ~ 5%. Their sources are the detector 
calibration, the aerosol component of the atmo- 
sphere and the statistical uncertainty of a fit to 
the HiRes/MIA data. A ~ 10% uncertainty in 
the fluorescence yield is common to both HiRes 
and HiRes/MIA and was therefore not included. 
The difference in the predictions of pure iron and 
pure proton maximum shower depths (X max ) be- 
tween different hadronic interaction models was 
not taken into account either, since we are not 
concerned about the fraction of real proton and 
iron showers here, but only about the fraction of 
showers with a certain X max . 

By generating two MC sets with pure proton 
and pure iron showers, we can calculate the effect 
a ±5% change in the proton fraction would have 



on the acceptance. The uncertainty in the final 
spectrum from such a variation is shown in Fig- 
ure[7] At the low energy end of the spectrum, the 
acceptance for iron cosmic rays is lower because 
iron showers develop higher up in the atmosphere 
and are more likely to be outside of our eleva- 
tion coverage than proton showers. This leads to 
larger uncertainties at lower energies. For ener- 
gies above 10 ls eV, where we use the HiRes stereo 
composition measurement, no difference was seen 
in the acceptances for iron and proton showers. 
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Figure 7. HiRes-2 energy spectrum with system- 
atic uncertainties (thick error bars) correspond- 
ing to a ±5% change in the proton fraction of the 
MC. 



5.3. Atmospheric Database 

Our current analysis uses a measurement of the 
average aerosol content of the atmosphere. We 
have studied the effect on the energy resolution 
of using an atmospheric database with hourly en- 
tries instead of the average. A MC set has been 
generated with use of the database to simulate 
data from September 2000 to March 2001. Our 
database covers 80% of all nights in this period. 
The MC set has been reconstructed first with the 
average and then with the database. A Gaussian 
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fit to the energy resolution has a a of 16.2% for 
the reconstruction using the database and 17.5% 
when using the average. 

An acceptance has been calculated in both 
cases. The ratio of acceptances can be seen in 
Figure [S] The acceptance for the MC set re- 
constructed with database is in the numerator, 
the acceptance for the same MC set, but recon- 
structed with the average, in the denominator. 
Using the database in the reconstruction does not 
have a significant impact on the calculated aper- 
ture. 
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Figure 8. Ratio of acceptances for a MC set gen- 
erated using an atmospheric database and recon- 
structed with the database (numerator) arid an 
average atmosphere (denominator). 



Data from September 2000 to March 2001 have 
been analyzed with the average and with the 
atmospheric database. We have not found any 
systematic shift in the reconstructed energies. 
When we examined the highest energy events (> 
10 195 eV), we found that using the database shifts 
one event from the bin centered at 10 1955 eU to 
the bin centered at 10 19A5 eV. This leads to a re- 
duction of the events above 10 19 5 from 7 to 6 in 
this dataset. 



6. Conclusions 

• Our measurement of the UHECR spectrum 
with the HiRes FADC detector in monoc- 
ular mode is consistent with the GZK flux 
suppression at 10 198 eU and with the Fly's 
Eye stereo spectrum. The "ankle" is ob- 
served at an energy of about 10 185 eU. 

• Detailed Monte Carlo simulation is used to 
calculate the acceptance of the detector, af- 
ter being tested with data-MC comparisons. 

• We avoid a bias of about ±20% in our spec- 
trum by using a realistic assumption about 
the input energy spectrum in the MC. 

• Systematic uncertainties caused by uncer- 
tainties of the input composition in the MC 
do not exceed our statistical errors. 

• We have tested the impact of using an 
atmospheric database with hourly entries 
rather than an average atmosphere in the 
reconstruction for about half of our ana- 
lyzed data. We have not found a significant 
difference. 
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